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In brane world scenarios with a bulk scalar field between two branes it is known that 4-dimensional 
Einstein gravity is restored at low energies on either brane. By using a gauge-invariant gravitational 
and scalar perturbation formalism we extend the theory of weak gravity in the brane world scenarios 
to higher energies, or shorter distances. We argue that weak gravity on either brane is indistinguish¬ 
able from 4-dimensional higher derivative gravity, provided that the inter-brane distance (radion) 
is stabilized, that the background bulk scalar field is changing near the branes and that the back¬ 
ground bulk geometry near the branes is warped. This argument holds for a general conformal 
transformation to a frame in which matter on the branes is minimally coupled to the metric. In 
particular, Newton’s constant and the coefficients of curvature-squared terms in the 4-dimensional 
effective action are determined up to an ambiguity of adding a Gauss-Bonnet topological term. In 
other words, we provide the brane-world realization of the so called i?^-model without utilizing a 
quantum theory. We discuss the appearance of composite spin-2 and spin-0 fields in addition to 
the graviton on the brane and point out a possibility that the spin-0 field may play the role of an 
effective inflaton to drive brane-world inflation. Finally, we conjecture that the sequence of higher 
derivative terms is an infinite series and, thus, indicates non-locality in the brane world scenarios. 


I. INTRODUCTION 

In the recent development of string/M theory |]p, branes have been playing many important roles. The idea that 
our universe is a brane in a higher dimensional spacetime has been attracting a great deal of interest . Although 
the idea of the brane world had arisen at a phenomenological level already in 1983 it is perhaps the discovery 
of the duality between M-theory and ifg x Es heterotic superstring theory by Horava and Witten [Q that made it 
more attractive. It actually gives the brane world idea a theoretical background: by compactifying six dimensions in 
the 11-dimensional theory, our 4-dimensional universe may be realized as a hypersurface in 5 dimensions at one of 
the fixed points of an SxjZi compactification. After cornpactification of 6 dimensions by a Calabi-Yau manifold, the 
5-dimensional effective theory can be obtained, see e.g. [||. 

Randall and Sundrum proposed two similar but distinct phenomenological brane world scenarios |^,||. In the 
scenario the 5-dimensional spacetime is compactified on SxjZi and all matter fields are assumed to be confined on 
branes at fixed points of the S 1 IZ 2 so that the bulk, or the spacetime region between two fixed points, is described by 
pure Einstein gravity with a negative cosmological constant. In the second scenario the fifth dimension is infinite but 
still the Z 2 symmetry is imposed. In both scenarios the existence of the branes and the bulk cosmological constant 
makes the bul k g eometry curved, or warped. There are generalizations of their scenarios with a scalar field between 
two branes In these generalized warped brane-world scenarios, the scalar field was introduced to stabilize a 

modulus called the radion, which represents a separation between the two branes. 

In the brane world scenarios with or without a scalar field in a warped bulk geometry, weak gravity in a static 
background has been extensively investigated Ml- It is known that weak gravity in the scenario without a bulk 
scalar between two branes is not Einstein but a Brans-Dicke theory at low energies. On the other hand, in scenarios 
with a bulk scalar field between two branes 4-dimensional Einstein gravity is restored at low energies. It is believed 
that the validity of Einstein gravity breaks down at a certain energy scale which can be much lower than 4-dimensional 
Planck energy. 

Hence, it seems natural to ask ’what does gravity in brane worlds look like at high energies or at short distances?’ 
In other words, ’how does the 4-dimensional description break down?’ 

In this paper we investigate weak gravity in brane world scenarios with a bulk scalar field between two branes at 
higher energies. Eor this purpose we use the gauge-invariant perturbation formalism developed in ref. |l^ . In this 
formalism all quantities and equations are Fourier transformed with respect to the 4-dimensional coordinates and 
classified into scalar, vector and tensor perturbations so that the problem is reduced to a set of purely 1-dimensional 
problems. We also adopt expansion in a parameter /x = where I is a characteristic length scale of the model 

and k^ is the 4-dimensional momentum (or the Fourier parameter). In the lowest order in /i, 4-dimensional Einstein 
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gravity is restored on either brane Pq . In the next order it is shown that gravity on either brane is indistinguishable 
from a higher derivative gravity whose action includes the Einstein term and curvature-squared terms. Equipped 
with the result for this order, we conjecture that in the order , gravity on either brane is indistinguishable from a 
higher derivative gravity whose action includes terms of up to the (A^-|- l)-th power of curvature tensors. Noting that 
the expansion in fj, is in principle an infinite series, this conjecture indicates that gravity on either brane is non-local 
at high energies even at the linearized level. This explains how the 4-dimensional description breaks down at high 
energies. Physically, the non-locality is due to gravitational and scalar waves in the bulk. 

This paper is organized as follows. In section H we summarize the basic equations by reviewing the formulation 
given in ref. [ p^ . In section III we perform the low energy expansion to investigate the system. In section IV we 
review linear perturbations in 4-dimensional higher derivative gravity so as to compare it with gravity in the brane 
world. In section ^ we discuss some physical implications. Finally, section VI is devoted to a summary of results. 


II. BASIC EQUATIONS 

The model and basic equations we shall investigate in this paper are exactly the same as those in ref. [p^ . We now 
summarize them briefly. 


A. Model description and background 

We consider a 5-dimensional spacetime A4 of the topology A44 0 IZ 2 , where represents 4-dimensional 
spacetime. We denote two timelike hypersurfaces corresponding to fixed points of the jZ^ compactification by E±. 
Each hypersurface can be considered as the world volume of a 3-brane. In order to describe we use the parametric 
equations 


Ei : x^ = (y±), (I) 

where {M = 0, •••,4) are 5-dimensional coordinates in A4 and each y± denotes four parameters {y^} (t* = 
0, • • •, 3). The four parameters play the role of 4-dimensional coordinates on each hypersurface. It is notable that the 
4-dimensional coordinates are not necessarily a part of the 5-dimensional coordinates. Actually, in the following 
we shall consider a 4-dimensional gauge transformation (a 4-gauge transformation) on each brane and a 5-dimensional 
gauge transformation (a 5-gauge transformation) in the bulk independently, where we call the 5-dimensional region 
between E± the bulk (and we shall denote it by Aib)- In particular, a quantity invariant under the latter (a 5-gauge 
invariant variable) is not necessarily invariant under the former (4-gauge invariant). 

We consider a theory described by the action 

Itot — I EH T ^matter : (2) 

where Ieh is the 5-dimensional Einstein-Hilbert action 

Ieh = 7^ [ d^Xy/^R, (3) 

^1^5 JM 


Iq, is the action of a scalar field 4' 


= — d^x\/—g 

JMb 


+ U{^) 


-Y. dSaV^Vai'fa) 




<7 — ± 


and Imatter is the action of matter fields confined on the branes 


dmatter — ^ ^ / d ; ^^tter]. 


( 4 ) 

( 5 ) 


Here, and q±fiu represent the pullback of and the induced metric on E±, and q±^,v is the physical metric on 
which is not necessarily equivalent to q±fj,v We assume that the physical metric is related to the induced metric 
by a conformal transformation depending on dtj.: 


q±fj,„ = exp[-a±(5'±)]y±^i,. 


( 6 ) 


2 








where a± is a function of 'I'j-, respectively. As shown in ref. |p^ the variational principle based on the action (|^) 
gives the correct set of equations of motion, including variations of , Qmn and dt. It is essential that the region of 
integration in the Einstein action (||) is not but Ad: the integration across gives the so called Gibons-Hawking 
term correctly. 

We consider general perturbations around a background with 4-dimensional Poincare symmetry: 


where the background is given by 


guN 

4' 

S±fJ,L' 


( 0 ) 

9 mn 
\|/(0) - 




(0)M 


- SgMN, 
54', 
-6Z^, 


S’il + 6S±^,, 


V]/(0) 

yi0)w 


c(0) 


2/±> 

w±, 

0 , 


(7) 


( 8 ) 


{x^} {g = 0, • • •, 3) represent first four of 5-dimensional coordinates {x^} {M = 0, • • •, 4) in Ad, w represents the 
fifth coordinate a:'^, and w± (w- < w+) are constants. Here, is the physical surface energy momentum tensor 
defined by 


= 


2 5 

^/-q± 5q±p.u 



d^ 2 /±>C± matter]. 


(9) 


and we have redefined V± and C± so that vanishes. Hereafter, we assume that the brane at ic = is our 
brane and that there is no excitation of matter on the other brane {S^'' = 0). 

The equations of motion for the background are as follows. 


and 


3A-|-3A^ = 
3A - 9i^ = 



( 10 ) 


( 11 ) 


where dots denote derivative with respect to w and 4*^^ = Here, we assumed that W- < w+ and that the 

bulk is the region w- < w < w^. The induced metric and the physical metric on for the background are 


„( 0 ) _ „-2A{w±) 

^± 11 V ~ ^ 

AO) _ 




( 12 ) 


B. Gauge-invariant variables 

Let us now construct gauge-invariant variables from the metric perturbation 5gMN, tbe scalar field perturbation 
54', the brane fluctuation 5Z^ and the matter perturbation 5S±^^. There are actually two types of gauge-invariant 
variables as there are two types of gauge-transformations: the 5-dimensional gauge transformation in the bulk (5-gauge 
transformation) 
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and the 4-dimensional gauge transformation on each brane (4-gauge transformation) 

( 14 ) 

As pointed out in ref. ||^ , these two kinds of gauge-transformations are independent. 

For the purpose of construction of gauge-invariant variables we expand all perturbations by harmonics in 4- 
dimensional Minkowski spacetime. This strategy is convenient since the background has 4-dimensional Poincare 
symmetry and the induced (and physical) metric on each brane is 4-dimensional Minkowski metric. In appendix]^ we 
define scalar harmonics T, vector harmonics V(t^l)ii and tensor harmonics T(^x,lt,ll,y)iii/- By using those harmonics 
we can expand all perturbations as 

SgMNdx^dx^ = {h(^T)T(T)^v + h(^LT)T{LT)fiu + + h(Y)T{Y)^v)dx^ dx'" 

h-ww Ydw^, 

= ipY, 

SZ±Mdx^ = {z±(T)V(T)^ + z±(^L)V(L)fi)dx^ -I- z±^Ydw, (15) 

and 

+ ^±(LT)T{LT)fj.u + ^±{LL)T{LL)iiv + T±{Y)T{Y)iivi ( 16 ) 

where we omitted dependence of harmonics and the corresponding coefficients on the 4-dimensional momentum 
and the integration with respect to The fc-dependent Fourier coefficients h(^T,LT,LL,Y)^ ^(T,L)tuj h^w, and "0 are 
functions of the fifth coordinate w only, and the other coefficients z±(^t,l )7 and t^(t,lt,ll,y) constants. 

Now we can analyze 5-gauge transformation of the coefficients of the harmonic expansion and construct 5-gauge- 
invariant variables, or those linear combinations of perturbations that are invariant under the 5-gauge transformation. 
The result is 


F{t) 

Fw 

F 

p 

■*- WW 

4>±{t) 

<t)±{L) 


h(T), 

h{T)w — e“^"^(e^"^/i(LT))’, 
h(Y) + ^AXu, + 

- 2e~^{e^X^y, 

z±{T) + , 

Z±{L) + ^(LL)L^^^ , 


4^:Lw 




Xu 


\w—W± ’ 


(17) 


where X^ = They form a maximal set of independent 5-gauge invariant variables constructed 

from the metric perturbation 6gMN, the scalar field perturbation 54' and the brane fluctuation SZ^. (The matter 
perturbation 5S±fiu is not included here since it is a 4-dimensional object.) 

We can also analyze 4-gauge transformation and construct 4-gauge-invariant variables, or those linear combinations 
of perturbations that are invariant under the 4-gauge transformation. For this purpose we first need to obtain 
expressions of various 4-dimensional quantities in terms of the 5-dimensional quantities gMN, 'k and Z^. What we 
need are the physical metric defined by the extrinsic curvature the pull back dtj. of the scalar field 

4', and the normal derivative 9 _l 4'± of 4> on I]±: 


'7±M>^(y±) 

K±tj.iy{y±) 

4'±(2/±) 

5j_4'± 


exp[ Q!±(4 'j-)] e^^e^j^gMNl ; 
^ M N r I 

2 ^±ti^±uFn±gMN\^^z±(v±) ’ 

'^\x=Z±{y±) > 


(18) 


where 


(19) 


M _ 


dy\ 


'± ' 


and £ represents a 5-dimensional Lie derivative. Using the harmonic expansions (|15D, we can obtain the corresponding 
harmonic expansion of 5q±fiu, 5K±^v, and The result is 

5q±^.v = ^±{T)T(T)fj.v + ^±(LT)T(LT)^iv + ^±{LL)T(LL)^iv + ^±{Y)T(Y)fj.i^, 

= k±(T)Ti^T))iu + k±I^LT)T{LT)iJ.u + k±(^LL}T(LL)tJ.v + k±(Y)T(Y)iiui 

= ip^Y, 

( 20 ) 


where 6K±^i, is defined by 




'(0)Pi 


( 21 ) 


and fc-dependent Fourier coefficients are 


v(«) 


o'±(T) = e ± F(t), 


_ ( 0 ) 

o’±(LT) = e “± </i±(T), 
_ ( 0 ) 

^±(LL) = e “± (?i±(L), 
(0) 

^±(y) = e 


F - 2A(I)±^, - kf,K(j)±{L) - a'±\e '^^ip + (()±„4'(‘’)) 


k±{T)=T-^e ^(e^^F’(T))', 

k±{LT) = ± 2 ®'^'^'^’ 
k±[LL) ^ 2 ^ 


1 


fc±(y) = ^2 |e-^(e^^F)- + (e^)'F^^ + 
1 


-e^r]^^''k^ki, - 2(e^)" 


4'± r, 


^/’±± = T 


2 L 


+ 2e^p + (t>±^ 


( 22 ) 

Here, the right hand sides of ( ^2[ ) are evaluated at re = w±, respectively, and have been written in terms of 5-gauge 
invariant variables Q The matter perturbation 6S±fj,i, on each brane can also be expanded by harmonics as 


SS±fj_L, = T±(T)T(T)^y + T±(lt)T(LT)pi> + T±(LL)T{LL)p,v + T'±(V)T'(v)pi/- 


(23) 


We can now analyze 4-gauge transformation of coefficients the harmonic expansion and construct the following 
4-gauge invariant variables from the physical metric perturbation 5q±^v 


JO) 


/±(r) = ct±(t) = e “± F(t), 

- _ 1 _ _ 

/± = CT±(y) + -'n^‘'kfj,k„a±(^LL) = e “± 


F - 2i(/.±^ - a'±^°^(e-2^^ + (/)±^4'(°)) 


(24) 


It is easily shown that k±(T ]^T,LL,Y)-: '0±i '4’±± and t^j’,lt,ll,y) are invariant under the 4-gauge transformation. 
Moreover, they are at the same time 5-gauge invariant since all except for are written in terms of 5-gauge 


^^The reason why they can be expressed in terms of 5-gange-invariant variables only is that they by themselves are 5-gauge- 
invariant 1^^. This fact illustrates that 4-gauge transformation is not a part of 5-gauge transformation and that these two 
kinds of gauge transformations are independent. 
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invariant variables and 5S±^i, is a 4-dimensional object. Hence, we have the set if±{T), f±, k±(^T,LT,LL,Y), V'±j 
^(T,LT,LL,Y)) of doubly-gauge invariant variables. 

From the point of view of observers on each brane, all observable quantities must be doubly-gauge invariant. 
However, all doubly-gauge invariant variables are not necessarily observable quantities. They can observe physical 
metric perturbation (/±(t)) /±) and matter perturbation only. The remaining doubly-gauge invariant 

variables (/E±(r,LT,LL,F)) '0 ±-l) shall be used to write down junction conditions of 5-dimensional quantities in a 

doubly-gauge invariant way. 

Our remaining task in this section is to give 5-gauge invariant equations in the bulk and doubly-gauge invariant 
junction conditions on each brane. Our final aim in this paper is to seek doubly-gauge invariant equations governing 
the physical metric perturbations and matter perturbations on our brane and to compare the resulting equations with 
the corresponding equations in 4-dimensional higher derivative gravity. 

Since there are many coefficients in the above harmonic expansions, let us divide these into three classes. The first 
class is the scalar perturbations and consists of coefficients of Y, T(^LL)^lv and The second is the vector 

perturbations and consists of coefficients of and T(^lt)ij.v The last is the tensor perturbations and consists of 
coefficients of T{T)^u- Perturbations in different classes are decoupled from each other at the linearized level. Hence, 
in the following we analyze perturbations in each class separately. Decomposition into scalar, vector and tensor modes 
is commonly used in cosmology. However, usually in cosmology we use scalar, vector and tensor representations of 
the isometry group related to the symmetry of 3-dimensional space. Meanwhile, here we will use scalar, vector and 
tensor representations of the isometry group of 4-dimensional space-time. 

In ref. it was shown that vector type perturbations vanish unless matter fields on the hidden brane are excited. 
In this paper we assume that there is no matter excitation on the hidden brane and, thus, we shall consider scalar 
and tensor type perturbations only. 


C. Scalar perturbations 


For scalar perturbations, we have three 5-gauge invariant variables from metric perturbation and scalar field per¬ 
turbation in the bulk: F, F^w and The first two are from metric perturbation and the last one is from scalar field 
perturbation. The Einstein equation leads to two relations among them 


and a wave equation 


F 

-i- II) 


-2E, 


2k24'(o) ’ 


A + 


(^ Y 

^2g3A 







= 0 , 


(25) 


(26) 


where a dot denotes derivative with respect to w and is the 4-dimensional momentum in the coordinate y^. 
Throughout this paper we consider modes with ^ 0 for scalar perturbations since a scalar mode with = 0 
preserves 4-dimensional Poincare symmetry and, thus, represents just a change of the background within the ansatz 
(^. Of course we shall consider modes with = 0 as long as 7 ^ 0. 

The boundary condition at ti; = w±, respectively, is given by the junction condition for the scalar held as 


T 




2e^if ■ 






(27) 


where t±(f) ^nd are doubly gauge invariant variables constructed from matter on and perturbation of 
the position of E±, respectively, and satisfy the following equations derived from the perturbed Israel’s junction 
condition 0 . 


^For a mode with = 0 we do not have the second equation of ( p^ since there is no tensor harmonics of the type {LL) for 
kfi, = 0. See appendix ^ for definition and properties of harmonics. 
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2 T±(y) = 

= TKie“'^““±’f±(LL). (28) 


Here, and is another doubly gauge invariant variable constructed 

from matter on I]±. The first of is nothing but the conservation equation of matter stress energy tensor S±f^i, on 
each brane. The second equation is the (LL)-component of the Israel’s junction condition and relates the perturbation 
of the brane position and the matter perturbation. The boundary condition (|^ ) at w = w±, respectively, can be 
rewritten to the following form by eliminating (p, t±(v) and 4>±w, and using the wave equation (|2^). 


C± 


F±2Kle 


—A—a 


( 0 ) 


(H + H")f±(z,z,) 






= 0 , 


(29) 


where 




w—w± 


(30) 


Finally, the doubly gauge invariant perturbation of the physical metric q±iiu is expressed as 


!± 


( 0 ) 

± 


F - 2A(l)±w 


„'^(0)(e-2A^ + ^i^vi,(0)) 


W — W± 


(31) 


In the next section we shall assume that 'i'(°)(zi;±)C'± ^ 0 to show the recovery of higher derivative gravity on a brane. 
In this case, the expression of f± can be rewritten to the following form by eliminating p and 4'±w, and using the 
boundary condition ([2^). 


JO) T 


e“S7± = F± 


(0) 


3a'- 

S7± I 






T±iLL) 


-W^ 


(32) 


Hereafter, we consider E+ as our brane and S_ as the hidden brane, and assume that there is no matter excitations 
on the hidden brane. Hence, we put = 0. 


D. Tensor perturbations 


For tensor perturbations we have only one 5-gauge invariant variable constructed from metric perturbation in the 
bulk: F^x) ■ The Einstein equation leads to the wave equation 


where k^ is the 4-dimensional momentum in the coordinate . 

The Israel junction condition leads to the following boundary condition at w = w±, respectively. 

[e^^F(T))' = ^±(T), 

where t±(t) is a doubly gauge invariant variable constructed from matter on Sj-. 

Finally, the doubly gauge invariant perturbation of the physical metric q±^i, is expressed as 

_ ( 0 ) 

f±(T) = e “± F(t) 


(33) 


(34) 


(35) 


Hereafter, since we assumed that there is no matter excitations on the hidden brane, we put = 0. 


7 















III. LOW ENERGY EXPANSION 


As already stated in the third-to-the-last paragraph of subsection [IB, our aim in this paper is to seek doubly- 
gauge invariant equations governing the physical metric perturbations and matter perturbations on our brane and to 
compare the resulting equations with the corresponding equations in 4-dimensional higher derivative gravity. Since 
we have only two gauge-invariant metric perturbations /+ and /+(r): and essentially two gauge-invariant matter 
perturbations 'f+(LL) and r+(T) on our brane ("f-i-cv) is related to t+(ll) by the conservation equation), we expect the 
following form of the equations on the brane. 


C{s)f+ = t+{ll), 

q^+^^''k^KC^T)f+(T) = T+(T), (36) 

where Cf^s) and C(^t) are functions of and is the inverse of the unperturbed physical 

metric The reason why we expect the linear dependence of metric perturbations on matter perturbations is 

that boundary conditions summarized in the previous section are linear in the matter perturbations. The reason why 
k^k^ was put in front of f+{T) is that we expect 4-dimensional gravitons on our brane (non-vanishing /+(t) 
with kf^ky = 0 and t+(t) = 0)- What is important here is that the functions C(^s,t) completely characterize the 
effective theory of weak gravity on our brane Since we are dealing with gauge-invariant variables only, there is no 
ambiguity of gauge freedom when we compare the result with the corresponding equations in 4-dimensional higher 
derivative gravity. Namely, we only have to compare functions C'(s,t) of k^^ki, with the corresponding functions 
of momentum squared in the Fourier transformed, linearized 4-dimensional higher derivative gravity. 

In this section we expand the basic equations summarized in the previous section by the parameter /i = k^k,y 
and solve them iteratively, where I is the characteristic length scale of the model which we shall determine by comparing 
the results of order 0(1) and 0(/i). The purpose of the ^-expansion is to analyze the behavior of the functions 0(5) 

and C{T) near /i = 0. Namely, we shall seek first few coefficients o|g ^) (i = 0,1, • • •) of the expansion 

OO 

q5.T)=E^”^(lT)- (37) 

2=0 


Since the 4-dimensional physical energy scale on S+ is given by m\ = —q^^^^'^k^kv = —fil 2 g“+W 2 A+^ 
expansion in /r is nothing but the low energy expansion. Hence, the first few coefficients o|g (i = 0,1 ,---) of 

the expansion (^) determine the low energy behavior of weak gravity on our brane. We expect that o|gj,j give 
4-dimensional Einstein gravity, that o|g give curvature-squared corrections to the Einstein gravity, and so on. The 

length scale I gives the energy scale l~^e°‘+ below which we can trust the ^-expansion. Nonetheless, we can 

defer the determination of I until we will obtain the results of order 0(/r) since I can be eliminated from all formal 
results in each order of the /r-expansion. Meanwhile, we shall keep it in intermediate calculations in order to make 
the expansion parameter dimensionless. 

First, by expanding F and as 

OO 

2=0 

OO 

nT){w)=Y.^i^Fl:^^{w), (38) 

2=0 

we can solve the wave equations (|^) and (^^ order by order. The result is 


^If one likes, one can restore gauge fixed equations for any gauge choices from the functions C(s t) only. For example, see 

( 0 - 
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( 39 ) 


F^^{w) = A{w)e^^^ 




dw 


A{w'y 


Jw- A(w'Y Jw- A(w") + A{w"Y 




pW 

of / + of 

J W — 


pw pw' 

+1-^ / dw'e^^^^''> / 

^ W— JW — 


with F^ where cW and Z3[*2 constants. 

Next, let us analyze the boundary condition. For this purpose, we expand f^(LL,T) as 


r+(LL) = 


2=0 


^+(r) 


+ (LL)> 


+ (T)> 


(40) 


i=0 


and put T_aAj) = ''’-(T) = 0 since we assumed that there is no matter excitation on the hidden brane. The boundary 
condition (E9[) for scalar perturbations can be rewritten as 


cf + B+cf = - 2^6“^^+”“+V'* 


+ (LL)’ 


where 


B. = 


B_ = 


x):' = 


cf + B_cf = 


r ,„4!w±dM,-3.,„, + fXl 

Iw- A^{w) Aj^ 


(41) 


A_ 

w+ 


d^ A(w)+A(wr e-3AM r 




Aiw)-^ 

g-3A+ ^ i(u;')e2^(™') 


J W- 

3e-^+ 


dw — 


A{w') + A{w')‘^ 


A(w') + A{w'Y 

f'^^{w') 

( 0 ) 


f[*1(w') 


\ 2 do) -Ai -aO’\4-,(0)-[*] 

r ^‘[w+) — K^a Ye + + w)|_ 


+ '+(LL) 


^ 3e ^ _ [,] 


(42) 


and X,£ = 0 . We have assumed that 4'^d± 7 ^ 0 and that A± ^ 0 . Here, in order to rewrite the last term of 
and x!)^, we have used the background equation We have adopted the following abbreviation: A± — H(w±), 

A± = A(r(;±), A± = A{w±) and = 4 /(o)(-u;^). xhe boundary condition (34) for tensor perturbations can be easily 
solved to give 

f ™ - 0 

= 0 , 


nW _ 

-^2 — rw+ 


1 _ [ 2,2 -2A+-a™-[*+!] _/-2w[*-l] 

r+ dwe-3AM r^(T) 

J W — 


(43) 
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where 


pWj. pw nw 

/ dw"e-^<-^'''^Fj!^^{w"), 

JW- jW- jW— 


( 44 ) 


and F[“^1 = 0. 

Thirdly, the doubly gauge invariant variables /+ and /+(t) corresponding to perturbation of the physical metric 
q+^iv on E+ are expanded as 


/+=£//W, 

2^0 

oo 

f+{T) =^d'ff(Ty 


i=0 


where the expansion coefficients are given by 




Sa' 


(0) 


2KlPa 


+ 




.(0) 


e“+ 4V) = Ct^(^+). 


iT)\ 


with fI 1] = = 0. 

Now let us summarize first few terms in the /i-expansion. 

/[°] = -167rGAr+n^f|°|^^), 


= -167rGAr+ 
f ™ - 0 

V(T) “ 




ft!rt = 


I+(T) 

/+(T) = 167rGjv+ 


r »2 , 2 -[ 2 ] ,,2 

+ ''T+'+(T) 




where 


= e 


167rGAr+ = 


+ . 


'■5“+ 


X 


= - 


JW — 

—3A(w) 


, . „ r+ , A{w) 

A{w) + A{w)‘^ Jw- A{w) + A{w)'^ 


+2 


^(■u;') + A(w')^ 


dw'e-^^^^"> 


3e-3^+ 

2^2^fG+ L 




1 + e^^+{2A^ ^ 


3e 


-3A_ 




^10+ ^10 ' /*^+ 

4+ X nA / dwe-^^^^'> = / dwe-^^^^'> / dw’e^^^'""> / dw"e-^^^^''\ 

Jw- Jw- Jw Jw — 


(45) 


(46) 


(47) 


(48) 


From these results, it is easy to show that 

1 + l's+q+^^''k^,K + 0(/i^) /+ = -167rGAr+f7+f+(ii), 

q+^^’"kf,K - lT+d +'’^''+ 0{A^)\ /+(T) = 167rGiv+T+(r), 


(49) 
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where q\_ is the inverse of the unperturbed physical metric, on S+. 

As promised, the length scale I does not appear in the results (49). However, we have obtained two lengths |^s+| and 
It+- It is evident that we cannot trust the /i-expansion, or the low energy expansion, at energies above min(|Zgj_|, 
Hence, we should choose I = max(|lg+l, /t+). 

The result (|4|) is of the expected form ( 3^ with the expansion (|3^) of the functions C'( s,t) ■ We have determined 
first two coefficients of the expansions: 


^[ 0 ] _ 

^(S) - - 

^[0] 


n 


-2 


/^[l] _ 

^(S) - 




167rGAr+ 

1 

167rGAr+ ’ 

= MO] 


y J 
It+ 

H-i-/ 


(S)> 


/^[O] 

U^T)- 


(50) 


IV. HIGHER DERIVATIVE GRAVITY IN FOUR-DIMENSIONS 

In this section, for the purpose of comparison, we review the linear perturbations in 4-dimensional higher derivative 
gravity. In particular, we consider a theory whose action includes curvature-squared terms: 

I = J [i? + , (51) 

where di, 02 and 0,3 are constants. The variation of this action plus matter action with respect to the metric gives 
the following equation of motion. 

G^i/ -|- diEifj^i, -|- d 2 E 2 ^i, + daAa^i, = SirGNTfii^, (52) 


where 


Eitj.v — -\- 2W^Rg^i, -I- 2RR^i, — 

E2fj.v = fj.^^Rup — iM^Rpp + W^R^y + Rpag^v + “^RppRy — 

Espy = 2{VPV^ + W^VP)Rppy„ + 2RpXpaRy^^'^ - \R"^aR''% 9 p.- (53) 

It is known that these three tensors Ei^y {i = 1,2,3) are not independent. Actually, it can be shown that Ei^y — 
^E 2 py + Espy = 0 in general. The easiest way to see this relation is to note that the choice (oi, 02 , 03 ) = (a, —4a, a) 
leads to a combination called Gauss-Bonnet term. Hence, the relation Ei^y — 4 E 2 py -I- Es^y = 0 follows from the fact 
that the Gauss-Bonnet term is topological. Because of the linear relation among Ei^y (i = 1, 2,3) it seems convenient 
to reparameterize the action di)- For later convenience we adopt the following reparameterization by taking another 
combination SR^R^ — R^. 

I = J d^x^g [E + aii?2 + a2(3i?0i?;) - R^) + 03 ( 1 ?'^^^%^^ " ’ (54) 

where ai = di — d 2 -I- 03 , 02 = 3 d 2 — dds and as = 0 , 3 . The equation of motion ( ^^ becomes 

Gpy + aiEipy + a2{iE2py — Ei^y) = SirGNTpy, (55) 

which is explicitly independent of 03 . 

We consider general perturbations around the Minkowski spacetime. Namely, we consider the metric of the form 

dsl = (g^°j + Sqpy)dyf^dyf^, (56) 
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where 


g(o) = (57) 

is a non-zero constant, and 

Sqfj.i, = ^{T)T{t)ij.u + ^{LT)T(LT)fj.u + ^{LL)T(^LL)fiu + ^{Y)T(Y)tJ.u- (58) 

Here, the coefficients cr( 7 ’), <^{y) are constants. For the definition of the harmonics T(^t,lt,ll,y)ij.vj see 

appendix^ Similarly to (|^), we can construct gauge-invariant variables /(t) and /. 

/(T) = d-(r), 

/ = d(y) -f ^ri^''kf,k^a(^LLy (59) 

As for the stress energy tensor we consider it as a first order quantity, and expand it as follows. 

Sij.1, = -I- f(LT)T(LT)nv + ^LL)T(LL)iii^ + T(Y)T(Y)fj.u, (60) 

where coefficients f(j’ ]^T,LL,Y) constants. In the Minkowski background these coefficients are gauge-invariant by 

themselves. 

We can expand the equation of motion up to the first order in the perturbations and express it in terms of the 
above gauge-invariant variables. The equation of motion is 


2r(y) = 3r]^’'k^k^T(^LL), 


1 -I- kfj,kv 


/ 




(61) 


for scalar perturbations. 


'^{LT) — 0 


(62) 


for vector perturbations, and 




kfj^kj. 


1 - 3a2q^°^>^''k^k, 


f(T) = IQt:GnT{t) 


(63) 


for tensor perturbations. 

Therefore, the weak gravity equations (|4^ (and the first of (^))on our brane is indistinguishable from the linearized 
gravitational equations and (1 ^ in the higher derivative gravity, provided that the following correspondence is 
understood. 


where Gn+j and are given by (^8|). 


Gn ^ Gn+, 

6 ai 

3u2 ^ 7 


(64) 


V. PHYSICAL IMPLICATIONS 


The result of the previous sections up to the order 0(/r) can be summarized by the following 4-dimensional effective 
gravitational action on our brane. 


/ = 


1 


167rGAr± 


J 


+ p 2 , M 2 


+ - Rl) + a±(i?f + 


-Rl) 


(65) 


where R±, R±i, and R^p^j are the Ricci scalar, Ricci tensor and Riemann tensor of the minimally coupled physical 
metric q±pi, on the brane S±, and constants Gat^, and are given by ( ^ ) (and corresponding expressions for 
quantities with the subscript “—”). The expression with the undetermined coefficient a± is the Gauss-Bonnet term 
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and does not contribute to the equations of motion at all. Hence 4-dimensional Einstein gravity is restored only at 
distances much longer than max(|^ 5 ±|, It±) or at energies much lower than min(|/gj_|, Otherwise, 4-dimensional 

Einstein gravity is not valid on the brane What governs weak gravity on the brane at higher energies is the higher 
derivative gravity. Note that the energy scale min(|^gj_|, can be lower than both 5-dimensional and 4-dimensional 
Planck energies. 

The restoration of the higher derivative gravity on a brane was almost insensitive to the form of potentials of the 
scalar field and conformal transformation to a frame in which matter on the branes is minimally coupled to the metric. 
Only condition which we had to impose is that there is a background solution with 4-dimensional Poincare symmetry, 
that 7 ^ 0 and that A± ^ 0 . The condition ^ 0 can be rewritten as follows ]T^ . 


V. 


'( 0 ) 




!l5y(0)y,(0) 


- -rVl 


(0)y^/(0) 


^0. 


( 66 ) 


This condition is equivalent to the condition that the background bulk scalar field is changing near the branes and 
that the inter-brane distance (radion) is stabilized. To see this it is enough to differentiate the background scalar-field 
matching condition (the second of (|ll|)) with respect to w±. In fact, the radion stabilization can be stated as 


0 7 ^ e 


^±. 


d 


dw± 




Xj>(0) 


±-14(4/^) 

w—Wj- A 


= + 






do) 


since the following equation is always satisfied because of the background equation. 


d 


dw± 




1 


w—w± D 


= 0 . 


(67) 


( 68 ) 


Hence, actually the condition (^) is equivalent to the condition that 7 ^ 0 and that the radion is stabilized. Since 
is related to 4'^°^ by ( 0 ) on branes, finally the condition ( | 6 ^ is equivalent to the condition that the background 
bulk scalar field is changing near the branes and that the radion is stabilized. The final condition A± 7 ^ 0 can be 
restated that the bulk geometry should be warped near the branes. 

A higher derivative theory with the action (|^) i s known to be unstable if the coefficient is negative ||^ . One 
can see this instability in the first equation of ( ph| ) at the linearized level: if Z|_|_ is negative then the equation has 
a tachyonic solution k^ki, = |Zs+|“^) with vanishing matter on the brane if^(LL) = 0)- Since the coefficient 

depends on the background solution as shown in (^^, one can perhaps conclude that the stability condition 
l's± ^ 0 should be imposed as a constraint on models of the brane-world. However, while we adopted the expansion 
in the parameter /r, the tachyonic solution corresponds to fj, of order unity. Hence, this solution is outside the 
domain of validity of the expansion in fj.. In other words, terms of more than the second power of /i can alter the 
stability/instability significantly. Further investigation is necessary to understand the stability of the brane world. 

One could also derive another stability condition at the linearized level from the second equation of (^). The 
stability condition at the linearized level would be < 0. Surprisingly, this condition is always violated since 

is positive as shown in (|^. Hence, the second equation of ( ^^ has a tachyonic solution k^^k^ = 

with vanishing matter on the brane (T+(r) = Oh The appearance of this tachyonic solution may be considered as 
the brane-world version of Horowitz instability 1261: noting that Horowitz instability is caused by the Weyl-squared 
term in a semiclassical effective action and that = (2/3)(3i?0i?)) — i?^) -I- (the Gauss-Bonnet term), 

one can expect that the term — R'^) in the action (|^) may cause an analogue of Horowitz instability. 

However, the tachyonic solution in this case is also outside the domain of validity of the expansion in /r. Hence, terms 
of more than the second power of /i can alter the instability significantly. From this point of view we expect that the 
tachyonic solution is just a spurious solution and does not indicate any physically harmful instability of the brane 
world. Actually, it can be shown that < 0 for any solutions of the full equation (^) insofar as the boundary 

condition (|^ with r±(r) = 0 is satisfied Q 2. The proof is easy: we can show that 


r'Wj^ 


rj^'^k^ki, I dwe^Fy,^ 


Wj^ 


dwe 


-3A 


[(e2^F(T))-] <0. 


(69) 


‘‘For a similar statement without branes, see ref. llC 
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Therefore, we conclude that the brane-world version of Horowitz instability is spurious and does not indicate any 
physically harmful instability of the brane world. 

Hence, we have to remove the spurious tachyonic solution for tensor perturbations whenever we deal with the low 
energy equation (^) or the low energy effective action (^). One of the possible ways to remove spurious solutions is 
the method of “self-consistent reduction of order”. This method has been used in many areas of physics including the 
radiation reaction equation the post-Newtonian equations of motion in classical relativity higher derivative 
gravity Q, and the semiclassical gravity |3^ (see also 

For some purposes it is useful to rewrite the gauge-invariant equations ( |4^ in a gauge-fixed form. We adopt a 
gauge (the harmonic gauge) such that = 0, where — dg±<5((/2 and 5q± = 6qf^^. Hereafter, indices are 


raised by 


= n 


-2 


do) 


and lowered by q±^i, 
is reduced to f± = (2/3)(t±(v), and the first equation of (^^ becomes 


= In this gauge the scalar-type gauge invariant variable f± 




1 -h I's^q^'^'' 




^±(y) = —1Q'!tGnt±{y)- 


(70) 


Hence, combining this equation with the second equation of i 
(P^) we obtain 


3) and the conservation equation (the first equation of 




+ l^±){d^diyO — — — 16 TrG]\[±Sf£^, 


(71) 


where is the surface energy momentum tensor of matter on the brane E± in the minimally coupled physical frame, 
□ = and h = = —Sq±. Alternatively, the equation (0 can be derived from the action (^) by using the 


well-known formula = —n[h^ — M^/2)/2 0 {h‘^) and ( p^ with ai = l‘s±l& and 02 = lx±/^- In the low energy 

limit the gauge-fixed equation of course, reduces to the linearized Einstein equation = —16t:Gn±S'^^. 

Again, we have to remove the spurious tachyonic solution for tensor perturbations whenever we deal with the low 
energy gauge-fixed equation & by eg. the method of “self-consistent reduction of order”. 

Given the high energy correction to 4-dimensional Einstein gravity, it is possible to derive corrections to the 
Newtonian potential. Actually, it is known that the spherically symmetric, static solutions of the linearized field 
equations in theories with curvature-squared terms are combinations of Newtonian and Yukawa potentials [l^ ]. Hence, 
the inclusion of the higher derivative terms indeed changes the short distance behavior of gravity, while it does not 
change the long distance behavior. For this reason, the correction presented in this paper is relevant for laboratory 
experiment on the validity of the Newtonian force. On the other hand, the long distance correction obtained by Garriga 
and Tanaka is relevant for astronomical tests although they considered a different model (the Randall-Sundrum 
infinite bulk model without a scalar field). Galculation of long distance corrections to the Newtonian potential in 
the model with a bulk scalar field between two branes is a worth while task as a future work. Analogy between 
calculations of long distance quantum corrections 0 and the brane world calculation of long distance corrections is 
also an interesting subject. 

Equipped with the result up to the order 0(/i), we conjecture that in the order 0{fj,^), weak gravity on the brane 
should be still indistinguishable from a higher derivative gravity whose action includes up to {N + l)-th power of 
curvature tensors. The coefficients of higher derivative terms can be in principle calculated by using the iterative 


results in section HI, as we have done in this paper up to the order 0{fi). Since the expansion in /i is in principle an 
infinite series, gravity in the brane world becomes non-local at high energies even at the linearized level. Actually, in 
this case the left hand side of ( 0 ) will have up to 2{N + l)-th derivatives {N oo) and, thus, we need to specify the 
0-th to {2N + l)-th time-derivatives {N oo) of the field hf^i, on a spacelike 3-surface in order to predict the future 
evolution. In other words, we need to specify the whole (past) history of to predict its future evolution, provided 
that can be Taylor expanded with respect to the time. This explains how the 4-dimensional local description 
breaks down at high energies. Of course, at low energies below min(|Z 5 ±|“^, l^].) we expect that the non-local behavior 
is suppressed. 

It is known that if a gravitational action includes a metric tensor (without derivatives) and its Ricci tensor only, 
then the system is equivalent to another system described by Einstein gravity and additional massive spin-2 and spin-0 
fields iQ. Evidently, this observation is applicable to the effective action (|65|). Hence, we can expect appearance 
of effective spin-2 and spin-0 fields on the branes. The appearance of these fields is consistent with the observation 
that the Yukawa potential appears in the correction to Newtonian potential by curvature-squared terms as a result of 
exchanges of effective massive spin-2 and spin-0 fields ||l^. In the brane world context, these spin-2 and spin-0 fields 
should be composite fields due to superpositions of bulk gravitational and scalar fields, respectively. 

The appearance of the composite spin-0 field may suggest a possibility of brane-world inflation without an inflaton 
on the brane. Actually, in 4-dimensional higher derivative gravitational theories it is known that inflation can be 
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driven by higher derivative terms. Such an inflationary model is called Starobinsky model [p2|-p4|. Hence, there is a 
possibility that the brane-world version of Starobinsky model may be caused by the higher derivative terms found in 
this paper. In this case, the composite spin-0 held is the brane-world version of the so called scalaron. This possibility 
of brane-world inflation due to higher derivative terms seems to be related to an interesting brane inflation scenario 
proposed by Himemoto and Sasaki pl|| . It seems interesting to explore relations between their model and the brane- 
world version of Starobinsky model. It may be sufficient to suppose that brane inflation is driven by the composite 
spin-0 field, or the brane-world version of the scalaron. It is expected that fluctuations of the spin-2 composite field 
as well as the spin-0 field (scalaron) play important roles in the generation of cosmological perturbations in this case. 

Although the composite spin-2 and spin-0 fields are generally expected to appear on the brane, it is not the end 
of the story. As far as the author knows, there is no known equivalence between a higher derivative gravity theory 
and Einstein theory with additional fields if the higher derivative terms in the action cannot be expressed in terms 
of Ricci tensor only. In other words, if the higher derivative terms depend on the Weyl tensor explicitly then the 
composite spin-2 and spin-0 fields are not enough to describe the whole system. Actually, the non-locality discussed 
above seems to require that the appearance of these composite fields is not the whole story since these fields have 
4-dimensional local actions. 

The expected non-locality should be due to gravitational and scalar waves in the bulk. Hence, the infinite series 
of higher derivative terms is one description of the non-locality pointed out in ref. ||3^ in the context of brane world 
cosmology 341 (see also |^). Another description was given as an integro-differential equation in ref. M, where 


a complete set of four equations governing scalar-type cosmological perturbations was derived by using the doubly 
gauge invariant formalism developed in refs. [ ^3p7| ]. One of the four equations is an integro-differential equation, 
which describes non-local effects due to gravitational waves propagating in the bulk. Further investigation of the 
relation between the two different descriptions may be an interesting future subject. 


VI. SUMMARY 

In summary we have extended the analysis of linearized gravity in brane world models with a bulk scalar field 
between two branes to higher energies by investigating the next relevant order in the expansion in the parameter 
^ where I is the characteristic length scale of the model given by I = 01,1^ max(|/ 5 ±|, Zt±)- Since 

the 4-dimensional physical energy scale m_|_ on our brane E_|_ is given by the expansion in /i is 

nothing but the low energy expansion. For the investigation we used the formalism developed in ref. [|l6[ , in which all 
quantities and equations including the surface energy momentum tensor and the junction condition are completely 
Fourier transformed with respect to the 4-dimensional coordinates so that the problem is essentially reduced to a set 
of 1-dimensional problems. 

We compared the result with the so called higher derivative gravity. It has been shown that in the order 0{p,), 
gravity on the brane is indistinguishable from the higher derivative gravity whose action includes the Einstein term 
and curvature-squared terms, provided that the inter-brane distance (radion) is stabilized, that the background scalar 
field is changing near the branes and that the background bulk geometry near the branes is warped. This result 
holds for a general conformal transformation to a frame in which matter on the branes is minimally coupled to the 
metric. The obtained indistinguishability between brane gravity and higher derivative gravity agrees with ref. [^8|] in 
which the non-relativistic limit of the theories with non-warped extra dimension was considered. In the present paper 
we considered weak gravity (including the non-relativistic and relativistic limits) in the brane world model with a 
scalar field in a warped bulk geometry. Newton’s constant and coefficients of curvature-squared terms except for the 
Gauss-Bonnet topological term have been determined as in (^). (The coefficient of the 4-dimensional Gauss-Bonnet 
term cannot and does not need to be determined as far as classical dynamics is concerned.) The result is summarized 
by the 4-dimensional effective action (|6^. In other words, we have provided the brane-world realization of the so 
called i?^-model. 

Equipped with the result up to the order 0(/r), we conjectured that in the order weak gravity on the brane 

is still indistinguishable from a higher derivative gravity whose action includes up to {N + l)-th power of curvature 
tensors. We discussed the appearance of composite spin-2 and spin-0 fields in addition to the graviton on the brane 
and pointed out a possibility that the spin-0 field may play the role of an effective inflaton to drive brane-world 
inflation. We also showed that the brane-world version of Horowitz instability is spurious. Finally, we conjectured 
that the sequence of higher derivative terms is an infinite series and, thus, indicates non-locality in the brane world 
scenarios. 
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APPENDIX A: HARMONICS IN MINKOWSKI SPACETIME 

In this appendix we give dehnitions of scalar, vector and tensor harmonics in an n-dimensional Minkowski spacetime. 
Throughout this appendix, n-dimensional coordinates are (/r = 0,1, • • •, n — 1), is the Minkowski metric, and 
all indices are raised and lowered by the Minkowski metric and its inverse rj^'^. 


1. Scalar harmonics 


The scalar harmonics are given by 

Y = exp{—ikpX^), (AI) 

by which any function / can be expanded as 

f = J dk cY, (A2) 

where c is a constant depending on k. Hereafter, k and dk are abbreviations of {k^} {p = 0,1,---,n — I) and 
respectively. We omit k in most cases. 


2. Vector harmonics 


In general, any vector field can be decomposed as 

= V(T)p + dpf, (A3) 

where / is a function and V(^T)ti. is a transverse vector field: 

= 0. (A4) 

Thus, the vector held can be expanded by using the scalar harmonics Y and transverse vector harmonics y{T)^. 
as 

V)J. = J dk [c(t)V(t)/x + C(L)dp.Y] . (A5) 


Here, c^t) and C(i) are constants depending on k, and the transverse vector harmonics V(T)/i are given by 

V(T)At = Upe^p{-ikpxP), 

where the constant vector Up satisfies the following condition. 

k^Up = 0 


for k^^kp ^ 0, and 


k'^Up = 0 , 

T^Up = 0 


(A6) 

(A7) 


(AS) 


for non-vanishing kp satisfying k^kp = 0, where is an arbitrary constant timelike vector. For kp = 0, the constant 
vector does not need to satisfy any of the above conditions. For the special case k^kp = 0, the second condition in 
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( [^ ) can be imposed by redefinition of C(i). Actually this condition is necessary to eliminate redundancy. Note that 
the number of independent vectors satisfying the above condition is n — 1 for yf 0 and n — 2 for = Q and 
that these numbers are equal to the numbers of physical degrees of freedom for massive and massless spin-1 fields in 
n-dimensions, respectively. _ 

Because of the expansion ( [A5| ), it is convenient to define longitudinal vector harmonics by 

= d^Y = -ik^Y. (A9) 


3. Tensor harmonics 


In general, a symmetric second-rank tensor field can be decomposed as 
where / is a function, is a vector field and t(T)/ii/ is a transverse traceless symmetric tensor field: 


— 0 




(T)tii 


= 0 . 


(AlO) 


(All) 


Thus, the tensor field can be expanded by using the scalar harmonics Y, the vector harmonics V(t) and V(i), 
and transverse traceless tensor harmonics as 


i/ii/ ~ j + C(^LT){d^V(T)y + 5i/^(T)/i) 

+ C(LL)(9^V(i)j, -I- i9jyV(L)^) + C(Y)Y'q^v'\ . 


(A12) 


Here, C(t), C(lt)i C{ll), and C(y) are constants depending on k, and the transverse traceless tensor harmonics Tfr') 
are given by 


= Sf_,„e:>q){-ikpxP), 

where the constant symmetric second-rank tensor s^i, satisfies the following condition. 


for k^kp 0, and 


k'^Sp^ = 0, 

sC = 0 


k'^Sp^ = 0, 

sC = 0, 


(A13) 


(A14) 




= 0 


(A15) 


for non-vanishing k^ satisfying = 0, where is an arbitrary constant timelike vector. For k^ = 0, the constant 
tensor Sau does not need to satisfy any of the above conditions. For the special case k^k^ = 0, the last condition 


in (A15) can be imposed by redefinition of C(^lt), C(ll) and C(y). Actually this condition is necessary to eliminate 
redundancy. Note that the number of independent symmetric second-rank tensors satisfying the above conditions is 
(n -I- l)(n — 2)/2 for k^k^ y^ 0 and n{n — 3)/2 for k^k^ = 0 and that these numbers are equal to numbers of physical 
degrees of freedom for massi ve an d massless spin-2 fields in n-dimensions, respectively. 

Because of the expansion ( A12 ), it is convenient to define tensor harmonics Ti^lt)^ and T(y) by 

T(LT)iiv = d(iV{T)i^ + dvV(T)p, 

— %{TLpkjj Ui/kp^Y^ 

T(LL)fiu = dfj.V(^L)i^ + d,^y{L)n - ViJ.i^d'’V(^L)p 


= -2kpku -I- -k^kpTipi, ) F, 


n 




(A16) 
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